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We have developed a model describing the non-proportional response in scintillators based on 
non-thermalised carrier and phonon transport. We show that the thermalization of e-h distributions 
produced in scintillators immediately after photon absorption may take longer than the period over 
which the non-proportional signal forms. The carrier and LO-phonon distributions during this 
period remain non-degenerate at quasi-equilibrium temperatures far exceeding room temperature. 
We solve balance equations describing the energy exchange in a hot bipolar plasma of electrons/holes 
and phonons. Taking into account dynamic screening we calculate the ambipolar diffusion coefficient 
at all temperatures. The non-proportional light yields calculated for Nal are shown to be consistent 
with experimental data. We discuss the implications of a non-equilibrium model, comparing its 
predictions with a model based on the transport of thermalised carriers. Finally, evidence for non- 
equilibrium effects is suggested by the shape of non- proportionality curve and wide dispersion in 
data observed in K-dip spectroscopy near the threshold. A comparison of the predicted curves shows 
good agreement for deformation potential value in the range 7-8 eV. 



I. INTRODUCTION 

The non-proportional response of light output versus energy deposition in scintillators presents a serious problem 
preventing otherwise good gamma-ray spectrometers from achieving a spectral resolution limited only by statistics. 
During the last decade considerable efforts were made in both experiment and theoretical modelling, resulting in a 
view that the non-proportional response originates from the non-linear interactions of electrons and holes in a tiny 
excitation volume leading to a quenching of luminescence^. The non-proportional part of the light yield is believed to 
occur during a very short (~ picosecond scale) period which we will refer to below as the non-proportionality stage. 

The diffusion of electrons and holes during the first few picoseconds after their generation is believed to play 
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a decisive role in the formation of a non-proportional response in scintillators. In all previous works referenced inS, 
it is assumed that hot carriers thermalise on a time scale of less than 1 ps and this assumption is central to all 
models based on the concept of thermalised carrier transport. While it may be valid for a few scintillating crystals 
it cannot be taken for granted a priori. In fact, it is possible that the opposite scenario, that is of non-equilibrium 
(non-thermalised) carrier and phonon transport, might actually be more realistic. 

The main objective of this work is to develop a dynamic model of electrons and holes in a high density neutral 
bipolar plasma interacting with longitudinal optical and acoustical phonons, and to use such modelling to analyse the 
quenching of luminescence in scintillators. 

In Section II we discuss the electron-phonon interactions in polar semiconductors, both for acoustic and optical 
phonons, taking into account screening in a neutral bipolar plasma. We derive an expression for the ambipolar 
diffusion coefficient at effective carrier temperatures exceeding room temperature, assuming that at high temperatures, 
even in extrinsic scintillators with a high activator concentration, the dominant scattering mechanism is due to 
electron-phonon interactions rather than impurity scattering. In Section III we derive and solve the balance equations 
describing the energy exchange between carriers, phonons and the thermal bath. In Section IV we evaluate the 
non-proportionality of response for Nal, (the scintillating material, for which most of necessary parameters can be 
found in the literature) and show that a non-equilibrium model gives results which are consistent with experiment. 
Finally in Section V we analyse the implications of our analysis and summarise the significant differences between 
non-equilibrium and thermalised models. 

II. ELECTRON-PHONON INTERACTIONS IN POLAR SEMICONDUCTORS 

In this Section we summarise textbook results for electron-phonon interactions and present them in a form which 
is convenient to use for the analysis of the non-proportional response in scintillators. There are several factors which 
we must identify as they enter the general expressions. They are: i) screening in a bipolar neutral non-degenerate 
electron-hole (e-h) plasma and ii) in the high temperature range of interest, T e 3> hQ^o, both the 2k < qsz and 
2k > qBz limits, where k is the carrier momentum and qsz is the phonon quasi-momentum at the Brillouin zone 
boundary, must be considered. 
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A. Screening in a non-degenerate, neutral, bipolar electron-hole plasma 



Following the absorption of an X- or gamma-ray photon, the carrier density in the tiny excited volume at the end 
of the photoelectron track can exceed 3 n > 1.0 • 10 20 cm~ 3 . At such a high density, ignoring the dynamical screening 
of the interaction between carriers and phonons leads to a serious overestimation of the rates of thermalisation, as 
well as the rates of momentum relaxation of carriers. In order to take screening into account we will use the Lindhard 



dielectric response function e(q, u>) for a non-degenerate semiconductor, e(q, uS) = 1 
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P(q, lu), where P(q, ui) is 



the polarisation operator, So is the static dielectric constant. Calculating the polarisation operator we write down the 
expressions for the real e'(q,uj) and imaginary e"(q,uj) parts in the form 
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Here to* and m* h are effective masses for electron and hole respectively, to is free electron mass, T is the carrier 



is the plasma frequency calculated for 



temperature (assumed to be equal for electrons and holes), Sl„ = 

V £om 

to* h = to. Other notations are x = - — — for the variable, /? = — - <C 1 and erfi(x) is the imaginary error function. 
In the expression for e'(x, lu) we took the limiting value e'(x, 0) having neglected small terms of the order of (3 -c 1. 



B. Phonon-electron interactions 



The absorption of phonons by a high density non-equilibrium e-h plasma acquires a special significance in relation 
to non-proportionality in a scintillator. Strong phonon re-absorption may slow down carrier thermalisation. After the 
absorption of an energetic X or 7 photon, the e-h density in the cylindrical excitation volume along the photo-electron 
track is very high, reaching for example a value of 2.0 • 10 20 cm -3 in Nal. Such a high density survives for at least a 
few picoseconds, which is crucial for the formation of a non-proportional response 2 ^ 3 -. The electron scattering times 
with respect to phonon emission and absorption are much faster (typically on a femtoseconds time scale), while the 
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FIG. 1: Phonon absorption rates by electrons and holes in NaT: (a) n = 2.5 ■ 10 20 cm 3 , (b) n = 1.0 ■ 10 20 cm 



diffusion of carriers as well as phonon propagation can be considered as slow processes (typical time scale of ps or 
slower). We must therefore discuss the evolution of the distributions of interacting hot carriers and phonons in much 
greater detail than was done previously. 

The rate of phonon re-absorption by carriers is one of the important parameters of energy exchange between the 
two systems. For longitudinal optical phonons (LO-phonons) interacting with carriers through the Frolich Hamilto- 
nian, the phonon scattering rate Tq due to interaction with carriers becomes 

r2 p"M 



Tq — -2M„— 



(2) 



Here u) = LUq is the dispersion relation for a phonon of a given mode with quasi-momentum q, Ma is the matrix 

element describing Frolich interaction in ionic crystals, M a = LO , e is the electron charge, — = [ J , 

e*T e* \eoo e a J 

£oo, is optical dielectric constant and hil^o is the energy of a longitudinal optical phonon. Finally P"(q, ui) is the 
imaginary part of the polarisation operator. 

We consider Nal as an example. Taking Eq = 6.6, — 3.0, TiQlo =20meV and m* = 0.287m, m* h = 2.397m 
as in Ref 6 , we plot the graphs of T/Qlo as a function of phonon momentum in units of k = \^2mhTf[o/h. Fig. 1 
illustrates the results. It is seen by comparing the results in Fig. 1(a) and (b) that screening has a substantial effect, 
reducing the strength of interaction for higher e-h densities, even though the carrier density in Fig. 1(a) is 2.5 times 
higher than in Fig.l.(b). 

Taking r/£lio=0.01 as a representative value, we obtain T v h-e = ~ 3ps. At the same time, the characteristic 
time for an energetic (~leV) electron or hole to emit a LO-phonon is estimated to be r e _ p /j ~5fs using the textbook 
formulas for unscreened electron-phonon interactions. This is 600 times smaller than r v h-e- However, electrons and 
holes are created with high excessive energies above the band gap edges. If we take the commonly accepted relation 
E e —h ~ 3E g , where E e -h is the mean energy for production of one e-h pair and E g is the band gap of the material, 
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then the mean electron (hole) excess energy is E g . This excess energy is rapidly converted into a population of 
LO-phonons, which cannot escape from the excitation volume because of their low group velocities. Even if the 
LO-phonon group velocity is taken as high as 1.0 • 10 5 cm/s, then before it is re-absorbed by electrons (holes) this 
phonon would have only passed a distance of ^3nm, less that the cross-section of the excited cylinder. This is one 
of the main reasons why LO-phonons do not escape from the excited volume during the whole period of the non- 
proportionality stage. This results in an effective bottle-necking of energy flow from hot carriers to the thermal bath, 
slowing down the temperature relaxation of carriers and justifying a model of non-proportionality in scintillators 
based on non-equilibrium transport. 

Another possible channel of relaxing carriers' excess energy is through the anharmonic decay of LO-phonons. 
The anharmonic rates of LO-phonons decay in Nal are difficult to estimate. In crystals such as GaAs LO-phonons 
are known to decay in about 5ps, although all decay pathways are open. It is expected that the anharmonic decay 
of LO-phonon in Nal will be slower. The reason is the large gap between the LO-energy and maximum energy of 
LA-phonons, %£Ilo > 2^0^^, prohibiting the possibility of LO-phonon decay into two acoustical phonons. For other 
scintillators, for example Csl, the same gap is smaller, thus allowing three-phonon anharmonic decay into acoustical 
phonons. However, most interaction channels will allow energy from LO-phonons to be transferred predominantly to 
TO-phonons, which have almost zero group velocities. Thus energy can escape from the excited volume only through 
the generation of acoustic phonons, being either the products of direct emission by thermalising electrons (holes) 
or of anharmonic decays of LO- and TO-phonons. In the first instance, the rate of energy loss due to emission of 
acoustical LA-phonons is typically one to two orders of magnitude smaller than that for LO-emission. In the second 
case, the energy released directly in the form of LA and TA phonons is only a small fraction of the LO-energy, and 
the TO-phonons contributing to the energy flow bottle-neck. 

Summarising, it is likely that over the duration of the non-proportionality stage the energy of the initial photon 
within the excited volume will be shared between hot carriers and LO-phonons. A larger fraction of this energy will 
be locked in the LO-phonon distribution. Indeed, in the fast process of LO-phonon emission, the excess energy of 
each carrier converts into the LO-phonon distribution. The excess energy released in the form of LO-phonons may 
be estimated as E g — T per carrier. For Nal E g =5.8eV, so taking T < 5000K implies that 5eV per carrier will be 
converted into approximately 2500 LO-phonons (20meV each). The actual numbers will be somewhat less if the 
parallel emission of acoustical phonons and possibly carrier-carrier interactions are taken into account. Even with 
1000 LO-phonons per carrier within the excited volume the following scenario is possible. Within a 5fs period each 
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carrier is likely to emit an LO-phonon; on the other hand the likelyhood of absorption of a single LO-phonon is 600 
times smaller (T p /i_ e =3ps). However, the probability of absorbing one of 1000 available LO-phonons per carrier is 
1000 times larger. This qualitative reasoning shows that over the non-proportionality stage there will be balance 
established between the carriers and LO-phonons. In a quasi-stationary state, the energy loss from hot electrons and 
holes will be balanced by phonon re-absorption. The hot carrier temperature, T, will be determined by the energy 
escape rate from the excitation volume through the emission of acoustic phonons and the diffusion of carriers leaving 
hot phonons behind. 



C. Scattering of electrons and holes by LO- and LA-phonons 



In this section we will analyse the rates of momentum relaxation of carriers due to Frolich interaction with LO- 
phonons, and deformation potential interaction with LA-phonons. General expressions can be found in textbooks 7 . 
We write down the standard expression neglecting dispersion for LO-phonons and assuming a linear dispersion relation 
for LA-phonons. We additionally take into account dynamic screening, and make no assumptions about the magnitude 
of the ratio 2k/qBz- 
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A e (h) is the deforma- 



tion potential constant for electrons (holes), p the crystal density and v s the longitudinal sound velocity. 

Taking into account only carrier interaction with LO- and LA-phonons we obtain the mobility of electrons (holes) 
(i e (h) m the form 
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FIG. 2: Schematic of energy exchange channels: (a) general scheme, (b) simplified scheme 



III. DYNAMICS OF HOT CARRIERS AND PHONONS 



In this section, we derive the balance equations for electrons and holes interacting with LO- and LA-phonons. 
Electrons and holes are assumed to be at the same quasi-equilibrium temperature at all times. This implies strong 
carrier-carrier (e-e, e-h and h-h) interactions. Therefore we may consider electron and holes as a single sub-system 
with a heat capacity equal to the total heat capacity of electrons and holes. The heat conductance describing the 
exchange of energy between carriers and phonons is a sum of individual (different) heat conductances, linking electrons 
and holes with the respective phonon sub-systems. Fig. 2 shows schematically the energy flow from the excitation 
volume into the thermal bath. Fig. 2(a) is more general, allowing energy relaxation of the system carriers+LO- 
phonons, where it was originally deposited, through direct emissions of LA-phonons and their subsequent escape 
(ballistic or diffusive) from the excited volume. Another relaxation channel is through the anharmonic decays of LO- 
phonons into their decay products, namely TO-, TA- and LA-phonons and subsequent escape of acoustical phonons 
from the excitation volume. Overall this balance scheme, although realistic, cannot be modelled accurately without a 
knowledge of the details of anharmonic interactions. For this reason we consider the simplified scheme depicted in Fig. 
2(b) where we incorporate all escape routes into LA-and LO-phonon decays, introducing an effective decay power 
PhA.decay{ThA) = 1laC L a{T L a - T ) and Plo, decay (T L o) = 1loClo(Tla - T ), where j L a(lo) are the effective 
decay rates, and Cla{lo) is the specific heat of LA- or LO-phonons. 

The balance equations describing energy exchange as depicted in Fig. 2(b) have the form 

dE 

= P c -Lo(Te,T LO ) + P c - LA {T e ,T LA ) + E c (0)S(t) 



dE LO 



dt 
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= -Pc-Lo(Tc,Tlo) - PLO, decay (Tlo, To) + 
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dE LA 



dt 



Pc-LA(T e , T L a) — PLA,decay{T L A, T ) (5) 



Here E c = E e + E^, E^q, E^a and T e , T^q, TlAi To are energy densities and temperatures of carriers, L0- 
phonons and LA-phonons respectively, To is the thermal bath temperature. In the right hand side of expression ([5]) 
Pc-Lo{Te,TLo) and P c -LA(T e , Tla) are the power densities dissipated by carriers into LO- and LA-phonons respec- 
tively, P c - L o(T e ,T L o) = P e -Lo(.T e ,T L o)+Ph~Lo(Te,T LO ) and P c _ LA (T e , T LA ) = P e - LA (T e ,T L A)+Ph-LA{T e ,T LA ). 
Finally, Plo, decay {TlOiPq) describes the dissipated power density which escapes from the excited volume via the 
anharmonic decay products of LO-phonons, and similarly PLA,decay(TjjA, Tq) has been introduced to describe the 
combined effect of LA-phonons being generated as a result of anharmonic decay of LO-phonons and LA-phonons and 
escaping from the excitation volume. Also in ([5]), E c (0) and -Elo(O) are energy densities deposited into electrons and 
holes and LO-phonons respectively, at t = 0, so that if Eq is the photon energy, then Eq = E c (0) + Elo(0), neglecting 
any direct emission of acoustic phonons in the process of the initial generation of electron-hole pairs. 



A. Power dissipation from carriers to phonons 



The general expression for the energy loss from the electron system due to interaction with phonons has the 
form 8 
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where is the electron distribution function, Nq is phonon distribution function, and is the electron dispersion 
relation. A similar expression can be written for holes. We will incorporate the additive hole contributions later 
directly in the balance equations. The expression for power dissipation in phonon collisions with the electron gas is 
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Here Mq is the appropriate matrix element for the electron-phonon interaction 
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The total energy of interacting carriers and phonons is conserved within the system leading to P e - p h(T e , T p h) + 
Pph~e(T e ,T p h) = 0. This conservation relation has already been used in the balance equations ([5|). 

Instead of evaluating the dissipated power P e ^ p h(T e , T p h) using ([6]) we write it as -P p h- e (T e , T p h) using expression 
([7]). For non-degenerate electrons this becomes 

P ph . e (T e> T ph ) = 4nh (J- - ±-\ cxp (Jf) x 
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where /i is the electron chemical potential. Calculating the integral over k we obtain 
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Substituting this result into ([9]) yields P e -Lo(^e, Tlo). Adding a similar hole contribution we obtain 
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Similarly for power dissipated into LA-phonons we obtain 
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where A is a mean deformation potential constant. 



B. Solving the balance equations 

The system of balance equations ([3]) has three unknown variables, namely T e (t),T^o{t) and TlaM (the three 

temperature model). It is convenient to introduce dimensionless variables, x(t) = T e (t) /T(0), y(t) = Ti,o(t)/T(0) 

and z(t) = TLA(t)/T(0), where T(0) is the initial temperature of carriers and LO-phonons, which is established after 

the generation of e-h pairs and quasi-equilibration within the carrier-LO-phonon system. We may also introduce a 

3tt 3 / 2 ph 4 

dimentionless time t defined in special units to = — =jt r , — . . The initial temperature T(0) can be found 
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from 

3nT(0) + N c T(0) = 2nE g (15) 

where N c is number of elementary cells per unit volume. This condition implies that the excess energy density, which 
will be released prior the recombination of the generated e-h pairs, is equal to locally thermalized (within carriers and 
LO-phonons) energy. At high temperature, 3nT(0) is the total energy of locally thermalised carriers in the excited 
volume, and N c T(0) is the energy density deposited into a single phonon mode (LO-phonon). 
Re-writing the balance equations in dimensionless units we obtain 
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Eq. ()16|) can be solved numerically provided that the main material parameters are known. We will test our model for 
three different scintillator crystals: Nal, Csl and ZnSe. For all these crystals most of important material parameters 
can be either found in the literature or estimated. Where we could not find relevant data we show parametric graphs, 
allowing the unknown parameter to vary within an expected range. 
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FIG. 3: Relaxation of temperature of carriers and phonons in Nal, n — 1.0 ■ 10 20 cm 3 : (a) - A = lOeV, "fLoto = "fLAto = 0.1; 
(b) - A = lOeV, 7lo£o = 0.1, 7la^o = 1; (c) - A = 5eV, 71,0*0 = iLAto = 0.1; (d) - A = 5eV, 7lo*o = 0.1, 7la^o = 1- Arrows 
indicate t=10ps 

C. Nal, Csl, ZnSe: solution of balance equations 



The estimate for the initial temperature T(0) in Nal from is T(0) = 4440K for n = 2.5 • 10 20 cm~ 3 . For 
n = 1.0 • 10 20 cm~ 3 we obtain T(0) = 2050K. In both cases, we have T(0) > T . Using p = 3.67g-cm~ 3 we obtain the 

"A - / ( T(0) ) S ' F ° r ^ = 10 e ^ we obtain to = 1.34ps for n — 2.5 • 10 20 cm -3 
and to = 1.97ps for n = 1.0 • 10 20 cm~ 3 . Figs. 3 and 4 show the solution of balance equations in Nal for the two 
values of on-axis e-h density, different values of phonon decay rates and two values of deformation potential constant 
(taking the same value for electrons and holes) falling into typical value ranges for known semiconductors, A = 5eV 
and A = lOeV. It is seen from Figs. 3-4 that within a lOps interval the decreasing carrier and LO-phonon temperature 
depends, as expected, on the effectiveness of LA-phonon emission through the value of A. For n = 1.0 • 10 20 cm~ 3 this 
means that the temperature decreases by ~ 50% for A = lOeV and ~ 20% for A = 5eV. For n — 2.5 • 10 20 cm -3 the 
decrease in temperature over lOps interval is larger, 60% and 40% respectively, mostly because the initial temperature 
is higher, and the LA-scattering dependence of temperature is stronger than for LO-scattering. 

Finally, it may look unrealistic to use the concept of a peak temperature of the distribution of LO-phonons as 
high as 4440K, exceeding the melting temperature (961K) of Nal . However, the distribution of LO-phonons created 
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FIG. 4: Relaxation of temperature of carriers and phonons in Nal, n = 2.5 ■ 10 20 cm 3 : (a) - A = lOeV, 7lo^o = 7la£o = 0.1; 
(b) - A = lOeV, fLoto = 0.1, jLAto = 1; (c) - A = 5eV, 7lo^o = ilaU = 0.1; (d) - A = 5eV, 7lo£o = 0.1, ^la^o = 1. Arrows 
indicate t=10ps 

by the energy relaxation of hot carriers in a tiny excited volume is transient and far from equilibrium, implying that 
lattice vibrations of only one mode are present, whereas all other modes remain undisturbed. The heat capacity of 
the LO-mode is only 1 /6 of the overall heat capacity, implying that the energy in this vibrational state or squared 
amplitude of vibrations would be l/6th of hypothetical equilibrium value at 4440K. In this way the amplitudes of 
atomic vibrations with only LO-mode excited to an effective temperature of 4440K are equivalent to the amplitudes 
of all vibrations at an equilibrium situation at 1/6 of 4440K or 740K, which is below the melting point. 

The estimate for the initial temperature, T(0), in Csl from (flip)) is T(0) = 2905K for n = 1.0 ■ 10 20 cm~ 3 . Using 
p — 4.51g-cm~ 3 and A = lOeV we obtain to = 2. Ops for n = 1.0 • 10 20 cm -3 . Fig. 5 illustrates the solution of the 
balance equations in Csl. The limits along the horizontal axis were chosen to cover the interval 0-10ps. In Figs. 5 (b) 
and (d) we modelled the effect of anharmonic decay of the LO-phonon, assuming the LO-phonon decay time to be 
lOps, falling within a typical range for the spontaneous decay of optical phonons of similar energy in semiconductors. 

For ZnSe, using the material parameters, we obtain T(0)=880K for n = 1.0 • 10 20 cm" 3 and T(0)=1700K for 
n = 2.5 • 10 20 cm -3 . These smaller values reflect the smaller band gap in ZnSe, E g = 2.82eV. The smaller value 
of the pair-production energy may well result in a higher on-axis density of e-h pairs within the cylindrical excited 
volume; thus the second estimate for e-h pairs density looks preferable. Nonetheless, we will also discuss the results 
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FIG. 5: Relaxation of temperature of carriers and phonons in Csl, n = 1.0 • 10 20 cm~ 3 : (a) - A = lOeV, 71,0*0 = 7la*o = 0.1; 
(b) - A = lOeV, jLoto = 0.2, j L At = 0.1; (c) - A = 5eV, 7^0*0 = jLAto = 0.1; (d) - A = 5eV, 71,0*0 = 0.85, j L At = 1 

for a smaller density, where the initial state is closer to equilibrium. Fig. 6 illustrates the solution of the balance 
equations in ZnSc for n = 2.5 • 10 20 cm -3 . The limits along the horizontal axis were chosen to cover the interval 
0-10ps (to = 3.1ps for A = lOeV). In Figs. 6(b) and (c) we use 7^Q=5ps, equal to that measured in GaAs, which has 
a similar phonon spectrum and close values for third order elastic constants. Also the data from literature* give A 
in the range 6.8-10.2eV. As seen from Fig. 6, even with a reasonably strong anharmonic decay over a lOps interval, 
the carrier temperature only falls to ^500K, remaining significantly higher that room temperature. For a smaller 
density n = 1.0 • 10 20 cm -3 without anharmonic decay the carrier temperature changes within <10% of T(0); with 
the estimated anharmonic decay it decreases down to 360K, implying that the carrier transport is not completely 
thermalised. 

IV. NON-PROPORTIONAL LIGHT YIELD IN Nal 

In this section we will discuss implications of the non-equilibrium transport model, checking its predictions versus 
the thermalised transport model of Li et aP . 

Following Payne et al we write down the scintillator yield in the form^£ 

Vpe = NoVcapVph (17) 
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FIG. 6: Relaxation of temperature of carriers and phonons in ZnSe, n — 2.5 ■ 10 2o cm s : (a) - A = lOeV, 7lo*o = Ila^o = 0.1; 
(b) - A = lOeV, -y LO to = 0.7, -jlaU = 0.1; (c) - A = 5eV, 71,0*0 = IlaU = 0.1; (d) - A = 5eV, 71,0*0 = 2.8, 7^*0 = 1. 



where Nq is the number of generated e-h pars, r\cAP is the efficiency by which the carriers energy is captured by the 
activators, and rjpH is the efficiency of light emission by activators. We will model t\cap using Birks and Onsager 
mechanisms^. 



r\cAP 



dE/dx 



(18) 



dE/dx\ B wks 



Here Tons is the Onsager radius defined according to 



= ksT, r e / h is the electron- hole separation, r/ e / h is the 



£0 r Ons 

fraction of free electrons and holes formed during the cascade, dE/dx is the electron energy deposition rate, and 
dE/dx\ Birks is an empirical fitting parameter relating to the strength of the exciton-exciton annihilation mechanism. 
The right hand side of (|18p is the scintillator light yield non-proportionality factor, which is a function of photon 
energy Eq and carrier temperature T, which we will designate as g(n). For a proportional scintillator g(n) = 1. 

Instead of the empirical Birks relation we will use a model of bimolecular quenching due to dipole-dipole Forster 
transfer and introduce a fraction of carriers surviving the non-linear quenching, 1 — QF, where QF is the quenching 
factor 2 . Thus, we obtain the local non-proportionality factor for a section of the track with carrier density n in the 
form 



g(n) = (1 - QF) 



1 - -q e/h exp 



ror. 



e/h 



gi(n)g 2 (n) 



(19) 



where Qi(n) = 1 — QF and £2 (ft-) = 



1 - rie/h exp - 
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To obtain the integral non-proportionality factor we 



r e /h 

must take an integral of (TT51) along the trajectory of photoelectron. If the quenching factor QF < 1/2 we may evaluate 
it analytically by integrating the quadratic form of linearised solutions to the two-dimensional diffusion equation with 
cylindrical symmetry of the excited volume and iterating up to the second order terms in the on-axis carrier density. 
As a result we obtain 

arctan (2y/ D a r/c 



£i(n) = 1 - ki^frn- 



kjrn 2 f 1 



2^/D^ 12 



arctan 



2^D a r/a 2 
{2^D a r/a- 



2 

+ 



arctanh - 



2./d r/a 2 arctann 
v " 1 ds Vv / 3T4l 



(1 + s 2 )\/3T4s 



2 



(20) 



where n is the on-axis density of e-h pairs, D a is the ambipolar diffusion coefficient in a bipolar plasma of electrons 
and holes, a is the radius of the excited cylinder, and r is the characteristic time, equal to the duration of an interval 
over which non-proportionality of response develops. The coefficient k2(t) — kij^ft is the bimolecular quenching rate 
parameter due to dipole-dipole Forster transfer. This has been measured in CsliS to be k2(t)y/t — 2.4- 10 _15 cm 3 s -1 / 2 . 
Following Li et aS^ we will use the same value to estimate exciton loss in Nal. For the ambipolar diffusion coefficient 
we take D a — = . For T e = Th we obtain D a = 



He + [l h Me + l-ih e fj, h 6 

Fig. 7 shows the results of calculation of the median temperature during the lOps interval as a function of e-h 
density. As seen from this figure the median carrier temperature during a lOps non-proportionality stage is much 
higher for A = 5eV in comparison with A = lOcV because of less efficient energy out-flow. 

Fig. 8 illustrates the dependence of the fraction of carriers surviving non-linear quenching as a function of density, 
calculated with the use of (J4j) and ([20]) . where we substituted values of carrier temperature for each density. The 
results shown in Fig. 8 were obtained by calculating the ambipolar diffusion from carrier mobilities using the Einstein 
relation. The mobilities were calculated accounting only for phonon scattering which is the dominant contribution 
at high temperatures. Any contribution from impurity scattering was neglected. With the increase in carrier density 
the fraction of surviving carriers Qi(n) decreases. With no other fitting parameters, except values of deformation 
potential, we obtain the curves shown in Fig. 8 corresponding to the middle of the range calculated by Li et al for 
fi e — fih within the range from to 800 cm 2 /V-s£. The closest to both curves shown in Fig. 8 is the curve from 
Li et al at /i=2cm 2 /V-s. However, the latter curve was calculated with the mobilities and diffusion coefficients of 
electrons and holes being constant throughout the whole range of density variation. For a non-thermalised e-h gas, 
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FIG. 7: Mean carrier temperature in Nal within lOps interval after photon absorption as a function of e-h track density: 
diamonds - A = 5eV, solid boxes - A = lOeV 
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FIG. 8: Qi(n) as a function of e-h track density: diamonds - A = 5eV, solid boxes - A = lOeV 



all transport parameters change independently. Fig. 9 illustrates the behaviour of the normalised ambipolar diffusion 
coefficient. It is seen that the ambipolar diffusion coefficient varies significantly. In addition, this variation is non- 
monotonic showing an increase in diffusivity in the high temperature region. A non-monotonic variation of the 
ambipolar diffusion coefficient occurs in the high density and temperature region and indeed these two factors may 
be of major importance for this region. High density results in stronger screening and increased carrier diffusivity. At 
high temperatures, the thermal momentum of carriers exceeds the maximum phonon momentum setting the upper 
limits of all integrals over the phonon distribution corresponding to the Brillouin zone boundary. This introduces some 
extra temperature dependencies into the expression for mobility An increase in the high temperature ambipolar 
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FIG. 9: Normalised ambipolar coefficient of non-thermalised carriers as a function of carrier density in the track, A = 5eV 

diffusion coefficients results in either saturation or an increase of Qi(n) over the high density region shown in Fig. 8. 

This feature is not present in any other models^— of non-proportional light yield and may therefore serve as an 

indicator of the non-equilibrated distribution of electrons, holes and phonons. 

Knowing the carrier temperature as a function of e-h track density allows us to evaluate the second factor in the 

expression for the activator capture efficiency, g2{n), describing the fraction of carriers surviving exciton formation 

while migrating towards an activator. The Onsager radius is a function of carrier temperature, which in general 

strongly differs from the bath temperature, as seen in Fig. 7. For the average initial electron-hole separation in the 

region of the track we choose r e / h = (2n) _1 / 3 . This is an assumption which is consistent with the currently used 

dE 

model of a cylindrical track of varying density. This expression for r e /h differs from r e /h ~ 3E g / —j— used by Payne 
et a££, which describes the mean distance that a photoelectron passes to deposit sufficient energy to generate a single 
e-h pair rather than e-h separation. Fig. 10 shows the results of our calculations. To plot this figure for Q2(E) as a 
function of photo-electron energy E instead of Q2(n) we used the conversion between E and the local e-h density in 
Nal derived by Vasiliev^. 

The shape of Q%(E) shown in Fig. 10 is solely connected with the non-equilibrium distribution of carriers through 
the dependence of the Onsager radius on carrier temperature. It is different in any other models of non-proportionality 
in scintillator response. In particular in the papers by Payne et a&£, rons is constant at T = 300-ftT, while r e / h 
decreases^ over the whole range 0.1-1000keV of E thus making the function 02(E) monotonic. The monotonic 
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FIG. 10: Q2{E) as a function of photoelectron energy, A = lOeV : dash - rj £ / h = 0.5, dash-dot - rj e / h = 0.8, dot - r) e /h = 1 



variation of the indepedent nonradiative fraction, "INF", was also reported by Li et a£. In their expression for the 
simulated local light yield the factor 1 — INF plays a role which is similar to Q2(E). Thus, both factors Qi(E) and 
Q2{E) have features which are specific only for the non-equilibrium scenario. 

Finally, combining all our results together we obtain g(E) as shown in Fig.ll. 

The "hump" structure is present in all simulations with different values of the parameter r) e /h- We see that as 
the exciton yield, 1 — ry e /fei increases the maximum value of g(E) decreases. However, the shape of the curves remain 
the same. The saturation of the non-proportionality factor g(E) at lower energies, or even the slight rise as seen for 
the case A = lOeV, is evident in both simulated curves in Fig.ll and experimental data below lOkeV. In simulated 
curves this is a specific feature of the non-equilibrium carriers arising from the special dependence of the ambipolar 
coefficient at high densities (temperatures) and corresponding features in g(E) at lower photo-electron energies. In 
fact particularly good agreement both for the shape and the amplitude of the curve on the left from its maximum is 
achieved in the range of deformation potential of 7-8 eV. The simulated curves deviate from experiment on the right 
of the maximum. We emphasize that we did not attempt the best fit, but simply took the important parameters a, 
t and &2 as they are quoted in the paper by Li et a£. For example, increasing a means that density of electronic 
excitations for any energy of photoelectron is smaller than what has been used in our modelling. This alone will 
result in a shift of the simulated curves in Fig.ll towards the left, because the same density will now correspond to a 
lower photoelectron energy. It is likely that the parameters ry e /^ and ki are functions of the non-equilibrium carrier 
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FIG. 11: q(E) as a function of photoelectron energy, T] e /h = 1: solid boxes - A = lOeV, diamonds - A = 5eV. Simulated 
curves are superimposed on experimental data of I.Khodyuk el (adapted from Fig. 11). Reprinted with permission from 
[I.V. Khodyuk, P. A. Rodnyi, and P. Dorenbos, J. Appl. Phys. 107, 113513 (2010)]. Copyright [2010], American Institute of 
Physics. 

temperature. In particular r) e /h (E) < 1 for higher E and carrier temperatures approaching room temperature. With 
density and temperature increasing, while E decreases, we expect r] e /f l (E) — > 1 for T 300K. Therefore, accounting 
for this dependence will result in an increase in the magnitude of g m ax as well as a corresponding increase in g(E) in 
the lower E part of the curves shown in Fig.ll. Overall, however, it will not affect the trend for saturation. 

Finally, the saturation or even rise of g(E) at lower energies is preceded by a change of curvature of g(E) on the 
lower energy side of its maximum. This change is not predicted by other models, although a change of curvature at 
lower energies was experimentally observed in Nal by Khodyuk et al 12 . Our modelling is consistent with the shape of 
their curve and the observed magnitude of non-proportionality, indicating that potentially the proof of non-equilibrium 
carrier effects was obtained in this experiment. 
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V. DISCUSSION 
A. Thermalised model 

To discuss the implications of a non-equilibrium model, we produce formal calculations of non-proportionality 
for the fast ( much shorter than lOps) thermalisation of carriers. We first check the degree of degeneracy of carriers 
thermalised down to room temperatures at high electron and hole densities. Even at n = 1 • 10 20 cm -3 , the electrons are 
strongly degenerate at To = 300K in all crystals. For all of them we have Ce/T) > 10, where £ e is the chemical potential 
of electrons. The Einstein relation for an arbitrary degree of degeneracy** is eD e ^ — [i e ^n — - — . Correspondingly, 

dCe(h) 

for strongly degenerate carriers the limiting value of the factor n is 2/3C, e (h), arL d the electron mobility will be 

an 

closer to the value that we previously estimated for temperatures exceeding 1000K, under the assumption of non- 
degeneracy. In any case the mobility of thermalised electrons remains much higher than that of thermalised holes 
with a moderate degree of degeneracy. As a consequence, the electrons have a small effect on the ambipolar diffusion 
coefficient, which will therefore be determined mostly by the holes. Heavy holes in Nal and Csl will have Ch — 0, so 

that the factor n-^- rs 1.4T making ambipolar diffusion faster than our estimate, based on the assumption of hole 
an 

non-degeneracy. Therefore, we will obtain the lower value of the non-proportionality of response if we use our results 
in expression (|2"0"|) taking T e = 300K. 

Fig. 12 shows the comparison of the factors Qi(n) for both thermalised and non-thermalised distributions. The 
upper curve was calculated for A-5eV. The curve for A=10eV essentially coinsides with it, and is not shown for clarity. 
The important point to note is that both curves show signs of saturation at high densities. However, the increase in 
Qi(n) for the high n end seen in Fig. 8 for the non-thermalised distribution at A=10eV is absent at room tempera- 
ture. This indicates that the increase in diffusivity at high track densities for A=10eV is due to extra temperature 
dependencies in the mobility of carriers arising from limitation of the integration over phonon momenta by Brillouin 
zone boundary rather than twice the thermal momentum of carriers. As seen by comparing the two curves in Figs. 12, 
if the thermalisation is fast and carriers thermalise before nonlinear quenching sets in, then the quenching factor 
QF = 1 — Qi(n) becomes small. This simply reflects the result of increasing the mobility of carriers when the temper- 
ature decreases. The second factor Q2(E) for this model is simply Q2{E) = 1 — ?7 e //j exp (—ro n s{To)/r{e/h(E)). With 
rons being fixed at room temperature Q2 (E) monotonically decreases with E as in other models 2 -*^ within approxi- 
mately the same range. To obtain a significant decrease in simulated light yield in this case the non-linear quenching 
must be pretty strong. This is not the case in our calculations. Fig. 13 shows the calculated non-proportionality factors 
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FIG. 12: Qi(n) as a function of track density, A = 5eV: solid boxes - thermalised distribution, diamonds - non-equilibrium 
distribution 
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FIG. 13: q(E) as a function of photoelectron energy, T] e /h = 1, A = 5eV: solid boxes - thermalised distribution, diamonds 
non-thermalised distribution 



within our model for both the thermalised and non-thermal cases. Overall the model of non-thermalised transport 
looks more realistic. 

The use of more definitive expressions for mobilities and diffusion coefficients, and the Einstein relation in general, 
surprisingly becomes important in view of recent work on ambipolar diffusion in the electron- hole plasma* 13 . It turns 
out that the ambipolar diffusion coefficient in a quasi-neutral e-h plasma differs from that derived from the Einstein 
relation in the situation when the electron gas cannot be viewed as ideal. The parameter which determines the role of 
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the interaction between carriers in a charged degenerate plasma is r s = (3/47rna 3 ) 1 / 3 , where as = h 2 £a/m*< h se 2 is the 
Bohr radius. With the densities of carriers in the photoelectron cylindrical track in the range close to 10 20 cm -3 we 
obtain r s > 1, and therefore the thermalised degenerate gas of carriers inside the track is non-ideal. For example, in 
Csl with m* = 0.235m and £o = 5.8 for n = 1.0T0 20 cm -3 we obtain r s ~ 1. Therefore, in a thermalised degenerate e-h 
plasma inside the track, the strong exchange and correlation results in a completely different ambipolar diffusion^. At 
much higher transient temperatures of a non-thermalised e-h plasma, the kinetic energy of non-degenerate electrons 
and holes exceeds the mean Coulomb energy and the e-h plasma will be close to ideal, validating the use of the 
standard Einstein relation. 

B. Temperature dependence of non-proportionality 

In recent experiments, Khodyuk et al 1 * demonstrated a strong temperature dependence of non-proportionality 
in LaBi'3:5%Ce for thermal bath temperatures ranging from 450 to 78K. Below we discuss this result using both 
thermalised and non-equilibrium models. 

At first sight the rapid thermalisation of carriers leading to an increase in their mobilities with decreasing 
temperature seems to be consistent with experiment. Thermalised carriers diffuse away from the excited volume 
faster, and non-linear quenching is less effective. However, thermalisation down to room temperature over a lOps 
time interval is already a problem unless crystals exhibit anomalously high anharmonic interactions. The problem is 
worse for thermalisation down to 78K. Below T ~ MIlo all dissipation rates become strong functions of temperature, 
and temperature relaxation slows down significantly. Another argument is that below room temperature the mobility 
becomes dependent on the level of impurities, saturating when the temperature decreases. With typical activator 
densities in extrinsic scintillators it would not be surprising to have mobilities saturating long before liquid nitrogen 
temperatures. Experimental data nonetheless demonstrates a strong variation in the non-proportional response in 
LaBr 3 :Ce from 300 to 78K(see Fig.l of refi^). 

A non-equilibrium model must also connect the observed dependence of the non-proportional response on ambient 
temperature with carrier mobility. However this connection is not direct, but comes through the potential dependence 
of varying temperature profiles in Figs. 3-6 on ambient temperature. The dependence on ambient temperature affects 
the relaxation profiles through the details of energy dissipation from LO- and LA- phonons in the excited volume. If 
the thermalisation rate is small, as in Figs.3(c)-(d), 4(c)-(d), 5(c) and 6(a) and (c), the profiles of relaxing temperature 
are not sensitive to the variation with ambient temperature. This is expected, because over a lOps time interval the 
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FIG. 14: Temperature evolution in Nal: n = 1 • 10 20 cm 3 , jlaU) = 1, A = lOeV. Upper curve - T = 300K, lower curve - 
To = 78K 

temperature of carriers remains well above T , so that changing T has a minor influence. However, in the situation of 
a more effective anharmonic decay, or higher group velocities for optical phonons allowing them to escape the excited 
volume, decreasing the ambient temperature means decreasing the relaxing temperature along the whole lOps interval. 
To illustrate how this may affect non-proportionality in Fig. 14 we show a comparison of the profile of relaxing the 
temperature for Nal with the profile calculated with the same decay and power dissipation parameters with the only 
the ambient temperature T decreased to 78K. 

As is seen in Fig. 14, the difference between T e (t) profiles for T = 300K and T = 78K reaches 200K by the end 
of a 5ps interval. Correspondingly, the effective diffusion coefficient for the lower curve is larger than for the upper 
curve. It is worth noting that the decay time j~[ A was chosen to be equal to t = 2.1ps. Ballistic LA-phonons in Nal 
pass a distance of 7nm on average over to. Thus, the effective "decay" rate of LA-phonons in the balance equations, 
jLAto — 1) describing the decrease of numbers of LA-phonons in the excited volume may be consistent with the flux 
of outgoing LA-phonons. The calculation of non-proportionality for the lower ambient temperature of To = 78K 
taking the effective temperature of carriers to be 15% lower than that of the scintillator at room temperature yields 
p(0.5keV)=0.98 as opposed to 0.89 for T = 300K - quite a significant variation. Clearly, the effect is dependent on 
material parameters describing energy exchange between the excited volume and its environment. 

The effect of bath temperature on the non-proportional response in doped alkali halide scintillators is of special 
interest to study experimentally. The exponential dependence on effective temperature in Q2 may significantly en- 
hance small variations in the non-equilibrium temperature of the carriers due to crystal cooling in comparison with 
scintillators, which do not exhibit the "hump" in non-proportional response, such as LaBr 3 :Ce observed by Khodyuk 
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et aU^. Under such conditions non-thermalised transport can cause an anomalously strong temperature dependence 
of non-proportional response in alkali halide scintillators. The experimental observation of such an anomaly could be 
a strong argument in favor of the non-equilibrium scenario. 

C. The effect of the non-equilibrium state on quenching characteristics 

So far we have discussed the implications of the carriers being non-thermalised during the non-proportionality 
stage by comparing carrier transport characteristics. Apart from transport properties, the effectiveness of non-linear 
interactions resulting in quenching of luminescence is determined by the dependence of various interaction constants 
on the mean energy of carriers and also on such parameters as rj e /h (discussed in section IV) determining the exciton 
yield. Therefore the concept of non-thermalised carriers in relation to non-proportionality in scintillators is more 
general than the concept of non-thermalised carrier transport. In particular, non-proportional response will depend 
both on the carrier diffusion coefficients and the non-linear interaction constants changing with time. 

It is well known that in equilibrium conditions the Auger recombination rates may show a significant dependence 
on temperature^, supporting our arguments. Similarly, the bimolecular quenching parameter /c2 (t) may also depend 
on varying carrier energy, which adds an extra transient variability to luminescence quenching. If such a dependence 
is present one should be able to experimentally distinguish between the scenario of distributions of non-thermalised 
carriers and that of carriers which relaxed their kinetic energy down to the temperature of the thermal bath, for both 
the same excessive non-equilibrated number. 

D. K-dip spectroscopy and non-thermalised electrons and holes 

The role of non-thermalised carriers in scintillator non-proportionality is most important in the parts of the 
photoelectron track with the highest density. In ordinary situations this is the part of the cylindrical photo-electron 
track closest to the end where ionisation is the most intensive. One promising new experimental technique for the 
identification of the non-thermalised carrier scenario, is K-dip spectroscopy^. In K-dip spectroscopy, the energy of 
incident X-rays, Ex, is chosen to be close to the ionisation energy of the K-shell, Ek ■ As a result the K-photoelectron 
is released with an energy of Ex — Ek, while the K-shell hole left behind may result either in emission of escape X-rays 
due to K a ^, L ai fj, M a ji transitions or Auger recombination, leaving holes in outer shells and generating the Auger 
electrons. Neglecting the escape peaks, we concentrate on energy deposition in the K-cascade, which is the process 
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of sequential Auger recombinations of the holes in K-, L- , M-... shells until holes start to be created in the valence 
band. In the K-cascade the first Auger electrons originating from K a< p, L a ^p and also M a ^ transitions possess an 
energy > 1 keV. They all are released at the absorption site but move in random directions giving rise to their own 
cylindrical ionisation tracks. Due to the randomness in the directions of the Auger electrons ejection, these tracks do 
not overlap. The highest density of the generated electron-hole plasma will be located at the end points of each track. 
The Auger electrons with energies exceeding 1 keV have low ionisation cross-section, so that the starting points of 
their cylindrical tracks are well separated from the absorption point. However, in the last stages of the K-cascade, 
the Auger recombination involves the outer N- and O-shells, and the Auger electrons are ejected with energies less 
than IkeV. These electrons will create a high density electron-hole plasma in the vicinity of absorption site. If the 
energy of the K-photoelectron is less than 1 keV, its e-h ionisation cloud will be also located in the immediate vicinity 
of the absorption site. For Ex 3> Ek the fraction of energy lost due to non-linear quenching in the vicinity of the 
absorption site is small and can be neglected. However in the K-dip spectroscopy experiment with Ex — Ek — > the 
highest e-h density region is created at the absorption site with the K-photoelectron ionisation cloud overlapping with 
those from the last few Auger electrons from the concluding stage of the K-cascade. t As a consequence, the density 
in this region may exceed the density of the e-h plasma at the ends of few cylindrical tracks. In this case the energy 
lost due to non-linear luminescence quenching in the central region may become comparable to, or even exceed, that 
lost at the ends of the cylindrical tracks. 

The experimental observation of an "end-point" , Ex — Ek — > 0, of non-proportionality in the K-dip spectroscopy 
experiment-is shows an increase in the dispersion of data points starting at Ex — Ek <lkeV. Khodyuk et al attribute 
this to the low accuracy of the measurement in this region, resulting primarily from the subtraction of two large, 
similar numbers and to systematic errors. We believe that this data spread is due to more fundamental reasons. 
Indeed, the degree of overlap of ionisation clouds in the vicinity of the excited atom (i.e., the absorption site) is 
subject to the randomness in the velocity distributions of the last few Auger electrons and the K-photoelectron. 
As a result, the actual density of e-h plasma in the vicinity of the absorption site will strongly fluctuate about its 
mean value giving rise to an additional spread in the data points - smaller non-proportionality for lesser overlap, and 
stronger for higher overlap. Subtracting the mean values therefore reveals the data spread due to these fluctuations. 

From an experimental standpoint, we strongly recommend that this spread is studied in greater detail. With the 
mean e-h density in the vicinity of absorption site substantially exceeding that at the end points of the cylindrical 
tracks, the expected peak temperature T(0) may show a significant dependence on the e-h density. Thus, the shape of 
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the end point (E x — Ek) of the non-proportionality curve in K-dip spectroscopy can serve to test for and potentially 
probe the effects of non-equilibrium carriers. Another possibility is to study the possible variation of this shape with 
ambient temperature, which will cause a change in the dynamics of hot carrier relaxation. 

VI. SUMMARY AND CONCLUSIONS 

We have developed a non-equilibrium model for the non-proportional response observed in scintillators. We 
complement previous models, based on the idea of fast carrier thermalisation and thermalised carrier transport, with 
a detailed analysis of energy exchange between carriers and phonons. The main conclusion drawn from this analysis 
is that unless the scintillator material possesses anomalously high anharmonicity, the energy flow from the excited 
volume containing c-h pairs generated by the absorbed X- or 7-ray photon, is impeded by the relatively immobile 
distribution of longitudinal optical phonons. Once the thermalisation processes have been studied in detail the best 
strategy for optimising scintillator response can be deduced. The present work indicates that the most promising 
direction to pursue will be to affect the thermalisation rate in order to create the conditions for faster out-diffusion 
of carriers from the initial high density ionization track. In practice, this may require either making thcrmalization 
faster or slower, depending on temperature dependence of carrier mobility and coefficients of Auger recombination or 
bimolecular dipolc-dipole Forster transfer mechanism of non-linear quenching. The rate of thermalisation may also be 
affected by creating a specific type of disorder and in this regard we note that scintillators need not be as structurally 
and stoichiometrically perfect as semiconductors to achieve good performance. 

We wish to thank I.Khodyuk for providing the data used in Fig.ll. 
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